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Abstract. We take an analytic approach towards the framework of parametric resonance and 
apply it on preheating and locked inflation. A two-scalar toy model is analytically solved for 
the Xcjy^x^ coupling for the homogenous modes. The effects of dynamic universe background 
and backreaction are taken into account. We show the average effect of parametric resonance 
to be that x's amplitude doubles for each cycle of 4). 

Our framework partly solves the broad resonance for preheating scenario, showing two 
distinct stages of preheating and making the parameters of preheating analytically calculable. 
It is demonstrated for slowroU inflation models, preheating is terminated, if by backreaction, 
typically in the 5th e-fold. Under our framework, a possible inhomogeneity amplification 
effect is also found during preheating, which both may pose strong constraints on some 
inflationary models and may amplify tiny existing inhomogeneities to the desired scale. For 
demonstration, we show it rules out the backreaction end of preheating of the quadratic 
slowroll inflation model with mass m ~ 10~®. For locked inflation, parametric resonance is 
found to be inhibited if (p has more than one real component. 
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1 Introduction 

Inflation was first proposed in [1] 30 years ago to solve multiple problems previously encoun- 
tered by the Hot Big Bang theory. It soon received public attention and subsequent works 
of others [2-9] have greatly perfected our understanding of inflation. Inflation has become 
so popular that almost all high energy theories have expressed their personal views on how 
inflation may be achieved, e.g. [10-14]. Even now, after 30 years of development, articles are 
still being entitled "inflation" on a daily basis. 

At the same time, reheating, the immediate subsequent scenario to inflation, was more or 
less in oblivion. It was not until 1994 that Kofman et al. came to realize[15] that the explosive 
particle production, which they named as preheating, should be first processed before the 
ordinary perturbative method is applied. Three years later, they proposed the analysis of 
reheating[16], covering numerous factors taken into account. (See [17] for a recent review.) 
Besides broad resonance, the narrow resonance also came into people's sight in the 1990s 
as a different non-perturbative particle production process during reheating[18, 19]. The 
analytic study of preheating however remains stalled afterwards, although numerical studies 
and applications on inflationary models have been going on anyhow. Despite that some 
applications of the theory in [16] have been quite successful, it is still very mathematically 
involved, and therefore, some of the underlying physics are not revealed. 

On the other hand, locked inflation[9] was proposed as a fastroll inflationary scenario. 
However, it was later found to be subject to multiple constraints, which together ruled out the 
whole parameter space[20]. Among them, parametric resonance is a crucial one which may 
terminate locked inflation much earlier than expected. Recently, locked inflation is recalled 
in a cyclic universe model which incorporated D-branes and curvature [21], bringing locked 
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inflation and locked inflationary contraction to each cycle. The proposed model however 
faces the same parametric resonance problem as locked inflation does. 

In this article, we construct an analytic framework of parametric resonance for the 
homogenous mode, in a way more straightforward and less involved. We then apply it on 
preheating and locked inflation, and discuss the possible inhomogeneity amplification effect. 
First, from section 2.1 to section 2.3, we solve a half cycle of parametric resonance and 
show the cause of the exponential boost effect. In section 2.4, we add up the effect cycle by 
cycle and give the exact solution to the Mathieu equation. We then put universe expansion 
into consideration in section 3.1 and investigate how backreaction acts as the terminator of 
parametric resonance in section 3.2. Such approaches then allow the discussion of preheating 
and locked inflation in section 4.1 and section 4.2 respectively. A brief discussion is carried 
out in section 5 on how parametric resonance may amplify the existing inhomogeneity and 
constrain inflationary models. The results are summarized in section 6. 

2 The Simple Framework 

To construct the framework, we consider the simple case first in a static background with 
Lagrangian density 

£ = -\m^ct>^ - \m\^ - Xct>\^ - \d^4>d^4> - \d^xd^X. (2.1) 

and X c'-re real scalar fields with bare masses m and M. Their interaction strength is 
characterized by the dimensionless parameter A. Here we consider the configuration of 
A|(/)p ^ IT? ^ Ax^ and ^ to neglect all backreactions. | | indicates the ampli- 
tude of oscillation when operating on a field. We allow to be negative, but still use 
to indicate the absolute value of for convenience when we compare it with other values. 
In many of our calculations however, we just neglect because it's too small. Throughout 
this paper, we will stick to the homogenous modes of both fields, so the spatial dependence 
will also be neglected in the notation. For simplicity, we will also use the word "parametric 
resonance" to indicate only the parametric resonance within our assumed parameter space. 
Please note this Lagrangian doesn't suffer from problems like negative infinite energy from 
the negative because it's just the effective part we take out from the total Lagrangian, 
so terms like the quartic ones don't distract us. As long as such terms are negligible, one 
doesn't need to worry about the validity of eq. (2.1). 

From the Lagrangian, we can write down the potential 

X) = \m^4>^ + \m\^ + M\\ (2.2) 

in which rr?^ jl dominates according to the above configuration. T^(</>, x) is visualized in 
figure 1(a), where the vertical lines also suggests the domination of m?(fP' /2. More details of 
the potential are displayed in figure 1(b) and figure 1(c). 

It is not difficult to imagine the motion of the two fields under such a configuration. 
Under <^ m?, (j)^s effective mass squared would be = m? + 2Xx^ ~ rn? and (j) will 
always be oscillating with frequency m and constant amplitude, unaffected by the motion of 
X- We can easily derive its equation of motion 

^ + (m2 + 2Ax')0 = O, (2.3) 
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(a) Landscape of V{<j), x) (b) Magnified with > 



(c) Magnified with < 



Figure 1. Contour-plotting of V{cj), x)- The first figure gives the overall landscape. The next two are 
zoomed in to the marked region and redrawn for details with different signs of M^. A darker color 
implies a larger potential energy. 



where dots mean derivatives w.r.t time. After neglecting x, we get the Simple Harmonic 
Oscillator (SHO) solution (p{t) = \cl)\ sinmt. 

x's effective mass squared = 2X<p'^ + ~ 2X(j)'^ is much larger than (p's most of the 
time, so it will oscillate much faster than <f>, as long as (j) is not very close to zero. When 4> 
crosses zero however, the motion of x would become very different from the stage of large 
(j). We will model the two stages in turn, trying to find out why and how x is boosted 
exponentially. Still, we write the equation of motion of x here for future reference. 



X + (M^ + 2A</.2)x = 0. (2.4) 



2.1 Rolling Stage 



The rolling stage is when cf) is not very close to zero, therefore x is oscillating much faster 
than 4>. It also requires that x's effective mass varies slowly and thus can be regarded as a 
constant in each cycle of x's oscillation. For most cases in which the rolling stage takes up 
more than half of (p's oscillation, one only needs to take into account the second condition, 
and then the first one is automatically satisfied. So the condition of this stage is 



dm?, 

A. 



<. 1. (2.5) 



Although this stage can be easily solved with simple approximations, we still would like 
to work on it again here, mainly to arrive at the variables we are interested in and confirm 
the rolling stage has a generally vanishing total effect on them. Suppose at time t = 0, x 
is at maximum with X = Ixl ^-nd X = 0; ^-^d scalar field (p is rolling down slowly with time 
derivative 4>. We take the effective mass constant and get the zeroth order SHO solution 
of eq. (2.4) 

X^°^ = Ixl cosm^^t. (2.6) 
We then derive the equation of motion up to first order x^^^ by introducing ^, 

X^^^ + mlx'^^^ + 4A#|x|t cos m^t = (2.7) 
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and get the solution with the initial condition x^^\^) = ^i^d x^^\0) = 0, 
At time t = vr/m^, when x^'^^ reaches maximum again, we have 



If we take when (/> is rolling down as an example, we will find out x gets an additional 
amplitude correction from the decrease of 0. To be precise, we consider the time interval 
At = Tr/m^, during which x rolls from one maximum to the other, x's amplitude is affected 
at a rate 



A|xl ^ 



1 m 



X 



At At ml 



\X\- (2.10) 



Replacing A with differential operator gives 



x|oc^, (2.11) 



which is the effect on |x| due to the variation of (/>. 

Similarly, we can calculate x's phase change resulting from the variation of (p, or 
directly. The motion of (p generates an additional x^^\7r/m^) which changes the time needed 
to push X back to its maximum again. The phase change is therefore 

^x = rt = — = 7;^oS—, 2.12) 

in which rriy^ and are the initial and final effective masses of x- 

The above calculation has shown how x is affected by (/>'s motion when (p is away from 
zero. Although the effects on x's amplitude (eq. (2.11)) and phase (eq. (2.12)) are strong, x 
can't get exponentially boosted in this way. This is because the effects during the upward 
rolling period of (p are exactly inverse from those during the downward period, so they 
generally cancel, x is increased when (p is rolling towards zero, but it get decreased when 
(p moves away from zero. The same thing happens to x's phase change, so the net effect 
vanishes. 

2.2 Zero-Crossing Stage 

When (p is small and (p is relatively large, the condition of rolling stage eq. (2.5) breaks, and 
the system enters the zero-crossing stage. The zero-crossing stage is the period (p crosses 
zero. During this stage, the motion of x is very much different from that during the rolling 
stage. It's because x's effective mass varies very fast in this stage and can become smaller 
than (/)'s effective mass when (p is crossing zero. 

Because ^ Ax^, we still adopt the approximation that (p is not affected by the 
motion of X) and ^ A|(/)p is still negligible. Here we set t = at (p = 0. In most 
applications of parametric resonance, \mt\ <^ 1 holds throughout the zero-crossing stage. We 
then get the motion of the two fields 

(p=\(p\mt, (2.13) 
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and thus 

X + 2X\^\^mH^X = 0- (2.14) 

Here we use r = t{2X\(j)\'^m?)i as the scaled time, and primes as derivatives w.r.t r. The 
equation of motion of x then becomes 

x" + r\ = 0. (2.15) 

The breaking of condition eq. (2.5) corresponds to |r| ~< 1. If one is only interested in 
the behavior of the system in |r| < 1, one can simply expand x Taylor series of r at r = 
and only take the first few terms. Numerical simulation tells us, however, the energy boost 
on X spans not only |t| < 1, but at an even broader range. To precisely calculate the boost 
effect, we need to take the analytic solution of eq. (2.15) 

xir)=oF^[l,-^^r')xo + oF^[l,~r')x'or, (2.16) 

where o-^i(a) z) is the confluent hypergeometric (limit) function, and the subscript denotes 
the value of variables at r = 0. 

Given the initial condition xo and x'o^ one can already calculate the effect of parametric 
resonance. In order to enable our further calculation, we take out the leading term of x for 
|r| ^ 1 (but \mt\ <^ 1 still holds). To do this, we expand oFi(a, —z) in series of 1/z. The 
leading term is 

oFi{a,-z)^ ^sm{-{2a + l)-2^/z)zi-2, ^ > 1. (2.17) 
\/tt V 4 / 



Here we want to calculate the energy boost effect from parametric resonance, so we 
characterize the boost rate with 

V^- Ijm log^^f^. 2.18 

The scaled energy density of x is defined as 

EsiT)^T^x\r)+x'\r), (2.19) 
with its relation with the (unsealed) total energy density of x 



E = X\<t.\Wx' + l^" = \l (2-20) 
Substituting eq. (2.16) and eq. (2.17) into eq. (2.18), we arrive at 



2 »V P(i)x?-2V2r(|)r(t)x„xi + 4r!(|h?. 

In the above equation, the limit r — t- +oo has already been taken. As a result, all 
non-leading terms of r are removed. The 2y/z in the sin function of qFi approximation is 
removed too. 

From eq. (2.21) we find the boost rate is hardly dependent on any parameter of the 
system, (m. A, etc,) or the energy density of any component, as long as our previous assump- 
tions stay valid. It is also easily verified by numerical simulations. This constant property 
is very helpful when we sum the boost rates across cycles of (f) in subsequent calculations. 
The fraction in the log function of eq. (2.21) should be of order unity, so the boost effect is 
significant and exponential. The boost rate rj acquires the same sign with Xox'oi which means 
the energy of x can either be boosted or dampened, depending on the sign of Xox'o- 
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2.3 Phase Delay 

In the above subsection, we have represented 77 with xo and x'oj the motion state of x at 
(/) = 0. The previous calculation however only gives half of the analytic solution — we can 
calculate rj with any given xo and x'o; but we haven't developed any method in deriving 
them two, except numerical ways. Therefore we will derive xo and Xq analytically in this 
subsection and provide a completely analytic solution. But before that, it's quite important 
to first clarify what "phase" is referred to here. 

The system has two degrees of freedom, (p and x, which means we need four independent 
parameters to fully describe the motion state of the system at any specific time. We are 
interested in the two parameters characterizing X) which can either be chosen as x and x'l 
or its amplitude and phase. From eq. (2.21), we have already learned the boost rate r] is 
a function of xo and Xo- Although x is not strictly a SHO at the zero-crossing stage, if 
we manage to represent the motion of x at </> = with an amplitude-like and a phase-like 
variable however, the amplitude-like variable would cancel and r] would then only depend 
on the phase-like variable. It's obvious rj may not stay constant in every half cycle of (/>, 
so reducing 77 to a single-parameter function allows us to track rj better when studying the 
boost effect across many cycles of 0. 

In order to find the phase representation of x at </> = 0, we attempt to solve xo and Xo 
reversely. It is reasonable to believe x is SHO-like at |r| ^ 1 and set the initial condition as 



x(r) = ^^^^sin/3, x'{r) = ^/Ej7)co^P (2.22) 

where /3 is the variable we use to represent the phase of x at r. We are only interested in (5 
here, so we come up with the combination 

|r|xcos/3-x'sin/3 = 0. (2.23) 

By putting in the solution eq. (2.16) and eq. (2.17), it simplies to 

r(5) sin(^ + y - = 2r(5) si„(^ + - I,),', (2.24) 

For a specific process with initial conditions given, xo and Xo should be definite and 
thus independent of the choice of r in eq. (2.24). For this reason, we contract the r^/2 into 
j3 and redefine j3 as 

P = /3(t) ^ /3o - + (2.25) 

where (3q is the variable for phase but it's now independent of time r. For short, we will 
however still use /3 to indicate /3o- In this way, eq. (2.24) also becomes independent of r, and 
simplifies to 

Xor(^) cos(/3 - J) = 2x'or(^) sin/3. (2.26) 
This immediately gives the solution 

XO = 2r(^)sin/3X, (2.27) 
X^) = rf?)cosf/3-^W, (2.28) 



4/ \ 4 
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Figure 2. xo and Xo functions of (5. x (in the label of y-axis) represents Xo and Xo the blue 
and purple curve respectively. The max goes over all j3 S [0, 27r) so all functions are normalized to be 
maximally at unity. The dashed curve is the expected Xo curve derived from xo if X were a SHO all 
the time. The purple curve's deviation from the dashed curve shows the phase delay in x's kinetic 
part at T = 0. Parameters are chosen as A = 1, to = 1, M = 10~'^, \(f)\ = 100, and they are chosen 
as such by default in the following figures. 

where X is a function of (3 and Es{t) and is not important in our calculation. 

By comparing it with the SHO solution, or the SHO-like solution eq. (2.22), one can 
find out the motion of x at r = is indeed non-SHO. The phase of its kinetic term is delayed 
by 7r/4. This is also confirmed numerically for all /3 G [0, 2it) in figure 2. Such a phase delay 
can be understood from the decreasing transfer efficiency from x's potential energy to kinetic 
energy as (j) is approaching zero. Given the (3 that makes xo = 0, i.e. x at the bottom of its 
potential at r = 0, such /3 doesn't induce the largest x'o any more. To understand it better, 
we compare such /3 with f3 = (3 + 5(3, where < (5/3 ^ 1. The difference caused by 5(^ is 
two-fold. On one hand, the phase /? makes x reach every value slightly earlier than phase (3, 
and such difference in timing gives (3 a larger effective mass of x when compare at the same 
value of Xj and in consequence, a larger energy transfer efficiency of order 6(3. On the other 
hand, for (3, x g^ts an extra up-climbing period just before r = 0. This effect decreases the 
energy transfer efficiency of (3, but only has order 6(3'^. As a result, the net effect is (3 has 
a larger energy transfer efficiency and therefore a larger kinetic energy. The maximum of 
kinetic energy is then shifted towards a larger (3, and its phase gets delayed relatively. 

We continue our calculation by substituting the set of solution eq. (2.27) and eq. (2.28) 
into the expression of r], eq. (2.21), and it becomes a very simple function of /3 



We are now able to compare our results with numerical ones. In figure 3, three curves of 
r] are drawn — the analytical result eq. (2.29), the numerical solution, and the semi-analytical 
result which is acquired by substituting the numerically derived xo and x'o into eq. (2.21). 
All three curves match perfectly. 




(2.29) 
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Figure 3. The energy boost effect w.r.t /3. All curves are horizontally shifted to be starting at 
rj = and dri/d/3 > when /? = 0, so they are aligned. The curves of analytical solution, numerical 
simulation and semi-analytic result overlap perfectly so only one curve is visible. 



So now we have successfully reduced the parametric resonance effect to a completely 
constant one — it depends only on f3, the phase of x- Ths boost rate is however independent 
of any of the parameters of the system, A, m, etc, or even the energy density of x- This is 
indeed the case from numerical simulations. Because the boost rate is independent of the 
energy of x, we can simply sum the boost rates when calculating the total boost rate for 
many cycles of <j). 

If we naively assume (3 is evenly distributed in [0, 2tt), and observe many enough inde- 
pendent half cycles of parametric resonance, the average boost rate can then be calculated, 
which is ^ 

(r/) = ^ ^ \(/3)d/3 = 0.346. (2.30) 

This actually means for every half cycle of 0, x's energy density becomes e^^^^ = 2 times of 
its original. For future reference, here we also define the maximum boost rate 

rim= max r/(^) = r/f|^) =log(l + ^/2) = 0.881. (2.31) 

/3e[0,27r) V8 / 

2.4 Long Time Evolution 

We now continue to study the long time evolution of x across cycles of 0. Suppose the initial 
phase of x is /?, by which we mean initially, when reaches its maximum, there is 

X = y^^^ sin /3, x = V2Ecos(3. (2.32) 

We can then use the boost rate r?(/3) acquired in the last subsection to evolve the system 
from one maximum of <j) to the next. To characterize the state of X) at every maximum of 
(p we use a polar coordinate system as the phase diagram of x- The angular coordinate is 
chosen to be /3, and the radial coordinate is e''", where r]a = ^ is the total boost rate. 

previous 
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Figure 4. Phase diagram of x at the neighboring maxima of 4>. The initial condition at the earher 
maximum of (p is set to a smooth unit circle with mark "-" . Sign "+" marks the later maximum, i.e. 
after half a cycle of Plots are colored according to their initial phase /?, so every "-" mark evolves 
to its corresponding "+" mark with the same color. The rotation of colors indicates the phase shift 
— the difference between the /3's at the neighboring maxima of (p. 

Therefore the x and y coordinates are proportional to x ^'^'^ X respectively. If we choose the 
initial condition as a constant energy density with evenly distributed /3, it then corresponds 
to a smooth circle with unit radius on the phase diagram. For a general overview, numerical 
results are given in figure 4, figure 5 and figure 6. 

These figures suggest a set of linear operations be applied on each "-" mark on the 
phase diagram, in order to get to its corresponding "+" mark with analytical methods. The 
linearity is also determined by the linear equation of motion of x- To clearly demonstrate the 
operations are linear, we represent plots on the phase diagram with (1,2) matrices {x,y)'^ 
of their Cartesian coordinates, and model the total evolution of field x between neighboring 
maxima of </) as a (2, 2) linear transformation matrix, which we call here as the transformation 
matrix. The transformation matrix should then be decomposed into two parts: 

1. During the zero-crossing stage, x is boosted with the boost rate rj a function of phase 
f3. So for an initial circle on the phase diagram, it is distorted to an ellipse because of 
the phase dependence of rj. For convenience, we rotate the ellipse by shifting (3 to align 
its major axis with x-axis. If the initial circle has a radius e^", we know at once the 
semi-major axis has length e''"^'''", and the length of semi-minor axis is e''""''"". The 
largest negative boost rate —rjm comes from the time-reversibility of x's equation of 



-9- 




+ 



+ 



+ 



+ 



+ 



+ 



+ 



Figure 5. Phase redistribution over Iralf a cycle of (p. This figure aims to show how the evolution 
over half a cycle of 4> affects the distribution of f3, i.e. the phase of x- It is based on figure 4, and 
derived by scaling every "+" plot to the same radius while keeping its phase /3 fixed. 

motion. When time-reversed, the largest (positive) boost rate ry^ always corresponds 
to the largest negative boost rate —rj^ and proves its existence. During this step, the 
plot in the phase diagram is transformed to (x',y')"^ which satisfies 



2. The ellipse is rotated with an angle which is generated from the difference of the periods 
of (j) and X- This is because the period of (j) is not an integer times of the period of x- 
Therefore x doesn't return to the exactly same position when <j) reaches maximum (or 
zero) again as where it was at (p^s last maximum (or zero). One can easily confirm the 
rotation angle is independent of phase /3, so the elliptical shape is preserved. Here we 
suppose the rotation angle that should be added to /3 is 6, so the plot is transformed to 



By definition, 6 is the phase x goes through from one maximum of (p to the next. So we can 
simply calculate the rotation angle 6 as 




(2.33) 




(2.34) 




(2.35) 




Figure 6. Redraw figure 4 and figure 5 witlr |0| = 99 instead of 100. 



in which the bare mass M has been neglected. Therefore the total transformation matrix is 



cos (5 — sin5\ / e''™ \ _ ( e^"^cos5 — e '''"sin 5 
sin 5 cos (5 / V ^ 6"'''" / ~ V ^^"^ ^i'^ ^ ^ 



where 6 = cosh r/^ cos (5 = ^/2cos5 and B is a (2,2) matrix of eigenvectors which is not 
important.^ 

We are mostly interested in the total boost rate during many cycles of (/>, which is 
represented by rja = log y^x^ + = log ^/ {x, y){x, y)'^ . Since (x, y)^ is multiplied on the 
left by the transformation matrix for every half cycle of (p, the total boost rate is actually 
determined by many powers of the transformation matrix. For this purpose, we have derived 
eq. (2.36) which allows us to only consider the power of the eigenvalues of the transformation 
matrix. 

For > 1 i.e. cos^ 6 > 1/ cosh^ r]m, we have b — Vb'^ — 1 < 1 < b + Vb"^ — 1. Because 
one eigenvalue is larger than the other, the overall effect is the system converges to the 
eigenstate with the larger eigenvalue. This means all initial phases converge to the same one 
after enough long time. Because the larger eigenvalue is greater than 1, x's total energy also 
increases gradually. 

For b^ < 1 i.e. cos^ 5 < 1/cosh^r/m, the two eigenvalues have unit magnitude and 
opposite arguments. Consequently there is no such phase convergence, and no overall change 
in x's total energy should be expected. 



'^To be more precise, one may want to calculate the transformation matrix via a three-step method — two 
rotations with angle 5/2 separated by a boost. The difference however only shows in the eigenvectors B, but 
the eigenvalues remain the same. 
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Such difference is also quite clear from the figures. The boost rate has dependence on 
the initial phase and transforms a circle in the phase diagram to an ellipse. Because of such a 
transformation, the phases converge towards x-axis. On the other hand, 6 rotates the ellipse 
and this may ruin the convergence if it's large enough. These two effects compete and, in 
consequence, the above two distinct behaviors of the system may occur. 

This also explains the exact solution of eq. (2.4), the so-called Mathieu equation. The 
solution of Mathieu equation has bands of stable solutions in which the amplitude of x 
generally remains at its initial scale as time passes, and bands of unstable solutions which 
have an exponentially growing |x| w.r.t time. Those bands correspond exactly to the non- 
convergent and convergent cases respectively. 

For better understanding, we can also calculate the phase that is converged to. The 
fixed points under transformation matrix eq. (2.36) satisfy y" /x" = y/x, i.e. 

y ^ xe^'''" sm6 + y cos 6 

X x" xe'^^i"^ cos 6 — y sin 6 

The solutions are 

y sinh rjm ± \/ sinh^ rjm — tan^ 5 



x e^™ tan 6 



(2.38) 



We can see how the solution is represented in the phase diagram in figure 7. From the 
figure we know at once we should take the minus sign in eq. (2.38) as the phase converged 
to. Based on these results, we can further derive x's long-time boost rate (which takes time 
far more than sufficient for convergence to take place, if there is). For the convergent case, 
the long time boost rate is simply the log of the larger eigenvalue in eq. (2.36), and for the 
non-convergent case, no boost occurs so the long time boost rate is zero. We thus define the 
long time boost rate as 



, cos^ S < K — - 



log ( cosh rjm cos S + y cosh^ rim cos^ 6 — l] , otherwise 



cosh rim 2 ^2.39) 



It's also visualized in figure 8. 

In figure 8, we can see the convergent-non-convergent boundary is at 7r/4, which is exact 
and can be confirmed from the expression of rim- This indicates only half of the 6 values give 
a long-time boost effect. If we consider the case with a 5 slowly scanning over [0,27r) at a 
constant speed, we would get the average long-time boost rate 



1 



2tt 



im) = TT- / m{S)d5 « 0.346. (2.40) 

^TT Jo 

Please note there is relation rii(—6) = ///(vr + 5) = r]i{6). 

It's quite interesting here to notice (r/;) "happens" to be very close to (r/). Although 
the expressions of and are very different, this however is not a coincidence. (77) and 
(rii) are actually identical and numerical integrals have confirmed that. It is not yet very 
clear why they are equal, but the following might be a reasonable explanation. Every value 
of 5 £ [0, 2tt) is actually hit the same number of times during the calculation of (77) and 
(rii). Consequently the transformation matrices of all 5's operate on the initial state the 
same number of times, although the order of the operations is different for (77) and {rii). The 



- 12 - 




Figure 7. The diagram showing the four fixed phases under transformation eq. (2.36), and four 
regions of other phases with different directions converge towards. Here 6 « 7r/6. The red line is 
the phase converged to and is thus fixed under the transformation. The dashed curves with arrows 
show the directions other phases converge towards in different regions. The blue line stands at 
the boundaries of the neighboring regions heading towards opposite directions, and therefore is also 
invariant under the transformation. The red and blue lines are eigenstates with the larger and smaller 
eigenvalues respectively. 

transformation matrices do not commute, but their commutation is only a rotation. The 
commutation could result in an either slightly larger or smaller boost rate, depending on 
whether the rotation is towards the major axis of the total transformation matrix or away 
from it. A huge number of such commutations is needed to make the two sequences the 
same. Some of the commutations contributes positively to the total boost rate while others 
negatively. When summed up, their effects cancel, and maybe one might be willing to think 
this as the reason of the equivalence between {rj) and {rji). 

3 Steps Towards Cosmological Applications 

So far, we have established the framework to consider the boost effect of parametric resonance 
in a static background. To apply it practically on cosmology, there are still a few steps should 
be taken first. One is the need to put our framework into the expanding background; the 
other is the consideration of backreaction. We will discuss them in this section. 
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Figure 8. The long time boost rate under different rotation angles S. Smaller 6 typically leads to 
larger boost rate because the eigenstate is more aligned with the phase of maximum boost rate. Large 
enough J's don't generate any long-time boost effect because they are the non-convergent cases. 



3.1 Dynamic Background 

When we apply this model in cosmology, a dynamic background is required.^ We will still 
neglect the backreaction in this section, so the effect of expansion on (p is simply a damping 
factor 1 01 oc a~2(^t), where a(t) is the scale factor. The equation of motion of x however 
becomes 

X + mx + 2X^\ = 0, (3.1) 
where H = da/adt is the Hubble rate of the universe, and can be regarded as a constant for 

3 

each half cycle of 0. Writing x = a^vj/ simplifies eq. (3.1) to 

^ + 2X(f)'^^ = 0, (3.2) 

in which we have neglected the effective mass contribution from universe expansion for 
A|(/)p ^ > H^. Therefore we can see the universe expansion also gives x an a~2 term. 

However the universe expansion has other impacts on our framework. As \(j)\ is damped 
by expansion, the rotation angle 6 no longer remains constant across cycles of (p. Therefore 
changes should be made to the calculations in section 2.4. For convenience, we use the 
notation to indicate the variable x of the nth half cycle of (p, or at the time when cp 
reaches its maximum the nth time. When all the subscripts (n) are equal in one equation, 
we may neglect some or all of them, as long as no confusion might be caused. Primes here 
are defined as the difference 

X(^) = X(n+i) (3.3) 

^Besides expansion, this model is also applicable in contracting phases without any prejudice. Because we 
are more concerned with the expanding phase, we will only discuss the effect of universe expansion here. The 
effect of contraction however, can be derived in the same way, and stages are simply reversely sequenced in 
most cases. 
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Different from the static version eq. (2.35), in an expanding background, (5(n) should be 
redefined as 

/■<(n + l)=i(n)+2^ 

(5(n) = / 2\/2A|(/>|(t)cosmMt (3.4) 

To finish the integral, we manually fix |(/>|(t) at 101(11)) so it gives the same result with eq. 
(2.35). This approximation however gives an error 



/■*(n+i) ^ 2V2A||</>l\nll 
A,5(,) < / 2\/2A||(/>|(„+i) - I I cos mMt = (3.5) 

•^*(n) 

where the outer I I is to take absolute values. For short, we have 



m 

Similarly we get 



AJ < (3.6) 



b' 


= ^(n+l) - S{n) = 


2\ 


/2A|0r 
m 


6" 




2\ 


/2AI0I" 






m 


AS' 


^ 2V2Xm"\_ 
m 


\6"\ 





(3.7) 
(3.8) 
(3.9) 



The universe expansion redshifts \(p\, giving 



lXl^ = e-^-l = , 3.10 

|(/>| 2m ' ^ ^ 

which consequently redshifts 6. The average effect of parametric resonance is then determined 
by the scales of 6, 5' and S" . 

When (5' <C 1, 5 changes very slowly compared with the oscillation of (p. The system 
therefore has enough time to converge to the eigenstate whenever there's any. Then the 
system can be regarded as trapped in its eigenstate which is slowly varying. The boost rate 
r]i{6) changes correspondingly as 6 slowly decreases, and the average boost rate in this case 
should then be taken as the long-time boost rate, (rji), defined in eq. (2.40). As S is slowly 
redshifted by universe expansion, the convergent and non-convergent cases are met in turn. 
This makes x's evolution stair-like at this stage. Half of the time, the amplitude of x enjoys 
flat platforms lasting for cycles of (p during which it doesn't grow at all. The platforms are 
separated by steep growing regions which fills the other half of the time. The two distinct 
stages come in turn to get the stair-like curve of x's amplitude. 

When (5' ~> 1 but 6" <C 1, S varies fast but 5' varies very slowly. In this case, 5 
changes significantly between neighboring half cycles of (j), but the amount changed stays 
almost fixed. For every half cycle of (j), an almost constant 5' ~ 1 adds to the initial rotation 
angle 5. Therefore for time long enough, the average effect is f3 has an even probability at 
each value in [0, 27r). This corresponds to an even distribution of /3 in r/(/3) and, results in 
the average boost rate which should be chosen as (r/) in eq. (2.30). 

When 6" ~> 1, the relevance between neighboring (5's or (5"s is small, so one usually 
considers as random in [0, 27r), being irrelevant with the history of evolution. When long 
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Figure 9. The numerical simulation result of many preheating processes compared with our statistical 
estimation. (See section 4.1.) Universal parameters are m = 10~^Mp, M = lO^^^Mp, and 7 = 10~^. 
The only varying parameter for different simulations is the interaction strength A, starting at 10"'* 
and ending at 4 x 10"'^, with a step increment 10%. The green plots are the data from numerical 
simulation, taken every time </> reaches maximum. The red line is the statistical average result and 
the two blue lines enclose the area within error lu, derived from A(?7(„)) in eq. (3.14). On the left 
half, we see a reasonable amount of plots lying outside the blue curves because we only expect 68% of 
the plots in la. On the right half the universe actually enters the small 5' stage and the error doesn't 
grow so much. However the blue lines are derived under the assumption that the whole evolution is 
the large 5' stage, so they are only functions of n and keep extending. Thus we find on the right half, 
the number of plots outside the blue lines is significantly decreased. The numerical data stop when 
backreactions start to be important. 



enough time is experienced, the random distribution of 6 gives the same results as the case 
of (5' ~> 1 and 5" <C 1 so {rj) should also be chosen as the average boost rate. 

In realistic cases however, it is rare to tell the last two stages apart. They share the 
same average boost rate {rj) and most of their characteristics are alike, so there is no need 
to distinguish them. Moreover, the second category has a 5" <C 1 but non-vanishing. To 
demonstrate the difference between the last two categories, the even probability of (3 needs 
to be shown. This however requires a prolonged constant value of 6\ which in most cases is 
very difficult to achieve with a non- vanishing 5" , even though it's much smaller than 1. For 
these reasons, we will treat them as a whole {5' ~> 1) in our following discussions and call 
it the large 6' stage, in which (5(n) is randomly distributed and independent of the history. In 
the same way, we name 5' <C 1 as the small 5' stage. 

We now fall back to check the approximation we made previously in eq. (3.6) and eq. 
(3.9), where |(/>| and are regarded as constants in every half cycle. They give errors 
A6 < \6'\ and A6' < \S"\, which fit perfectly with the three stages above. Therefore this 
approximation is applicable here. 

Before we end this discussion on dynamic background, we would like to briefly calculate 
the statistical results from the random 6 for the large 5' stage. In ^ cycles of (p, the boost 
effect of parametric resonance occurs n times. The total boost rate is 

n 

^(n)('J) = E^('^(i))- (3.11) 

1=1 
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Assuming absolute randomness of all E [0, 27r), we get the expectation value of ri[n){^) 

(^(n))^^ I d"<5fjr?(5(i)) = ^£\K5)d5 = n(r?)« 0.346 n. (3.12) 

This is obvious because it's just how {r]) is defined. To further compute the statistical error, 
we calculate 

= n(r/2) +n(n- l)(ry)2. (3.13) 
and thus the statistical error can be estimated as 

^(^(n)) ^ ^(^(n)) - (^(n))' = V^a/^T^ « 0.539^^. (3.14) 

To justify both the statistical calculation and the definition of the large and small b' stages, 
here we also give figure 9 through multiple numerical simulations. 

3.2 Backreaction 

In this section, we discuss the backreaction of x by taking into account the Ax^ term in m|, 
the effective mass squared of (/>. Since we assumed homogeneity, we will not discuss other 
backreactions here. In most cases of parametric resonance, the assumption Ax^ ^ starts 
valid. However as x's amplitude gets boosted by resonance, this assumption would finally 
break at |x| ~ m/\/A, provided the boost effect hasn't yet been halted by other causes. 
Because |xP oc ra~ 1 « l/V2A(/.2, for every half cycle of ^, |x| reaches its maximum when 
— 7- 0. From this point of view, the violation of this approximation first starts small — only 
at — and then expands to larger (/> gradually while |x| is being further boosted. 

From the relation |xP fx 1/y^2A^, we can find the effective potential for (/>, which is 
actually 

Kff(</') = \vr?^^ + A</)'x' = \n?<\>^ + \Ax?m\mW. (3-15) 

where |x|m is x's amplitude at maximum (/>, and | |abs means taking the absolute value 
(and I I still indicates taking the amplitude). Here we can see the additional part of T4ff((/>) 
is proportional to the absolute of (/>. Compared with ^rr?^ which is proportional to 0^, 
the effective potential from interaction 14ff((/>) agrees with our previous instinct that the 
backreaction first comes to important at small 0. 

So we consider the motion of for half a cycle, starting from t = 0, = —\(\)\^ </> = 0, 
and Ixl = Ixim- We will use 7^ = Ajxlm/"^^ simplicity. When < 0, (/)'s equation of 
motion is 

+ m2<^-A|x|™|</'| =0. (3.16) 
We can analytically solve this equation and arrive at = 0. The solution is 



t = ti = 



(3.17) 



^(t = h) = (7^ + l)m|^| sinmii. (3.18) 

At (/) = 0, X is boosted by parametric resonance. This gives the interaction part of 
Ves{4>) a factor of e'^^. Here r] is still the boost rate. For this reason, the equation of motion 
is slightly modified for > 0, 

^ + mV + A|x|^|0|e2'' = O. (3.19) 

We then solve it again to arrive at = and finish the half cycle. At the > side when (f) 
reaches maximum again at time t = t2, we have 

TT 1 ■7^6^'' 

^2 = ii + arctan / (3.20) 

2m m ^ I ^ ^ 

^{t = t2) = mVl^e^"" + 272 + 1 - 72e2^'), (3.21) 

and 

\x\=e^\x\J-^- (3.22) 

We can clearly see when r] = 0, there is no boost effect and the motion of the system becomes 
symmetric w.r.t = 0. 

When the backreaction is weak, i.e. 7^ <C 1, we can expand the results to first order 
and get 

(3.23) 

m m 

0(t2) = |0|(l-(e2^-l)7'), (3.24) 
Ixl =e''|xU(l + (e'''-l)7'), (3-25) 
6^'^{l-{lie^^ + l)-2y). (3.26) 

Therefore, in case the backreaction presents, the oscillation of (p becomes faster and that of 
X is compared slower. The non- vanishing 77 gives a boost effect. Amplitudes (/>(t2) and 1x1(^2) 
are thus affected and become slightly different with their original values at t = 0, suggesting 
an energy transfer. 

Whenever one takes this backreaction into account, one should adopt these corrections. 
However for weak backreactions in which 7^ <c 1, the corrections actually can be neglected. 
Such backreaction only comes to important when 7^ ~> 0.1. In some cases, however, the 
exponential boost effect is terminated by other causes before 7^ reaches 0.1, so it never comes 
to important. 

If the system does evolve to 7^ ~> 0.1, one then needs to take it seriously. The most 
significant effect is oscillates faster and faster, due to the effective mass contribution from 

the backreaction. This has several consequences. First and most obviously, because the 
boost rate is counted per half cycle of (p, a faster rate of 0's oscillation also leads to a faster 
boost to X- The faster boost acts back on 0's effective mass and causes to oscillate even 
faster. The positive feedback gives an explosive boost effect to x once the backreaction 
comes to important. The explosive boost effect is soon shut down when \x\m reaches the 
same scale with \(p\. This is because (p and x would then have similar frequencies and the 
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approximate replacement of with (x^) would be no longer applicable. This immediately 
ends the exponential boost effect of parametric resonance. 

Moreover, this backreaction may also offer a major contribution to 6' . Here we give 
a simple estimate based on eq. (3.26). By definition, we have (7^)' = 2A(|xlm)7"i2 = 
^2^g2?7 _ go the correction to 6' up to first order of 7^ is 

A6' = _2%^2W/vr 4,^,2,' _ _ ^2, + (3.27) 
m V 2 / 

When it has the same magnitude with the originally defined 6' in eq. (3.7) in an expanding 
universe, we can solve the magnitude of 7^, 

^2 < 0.1. (3.28) 

m 

This means the impact of backreaction on rotation angle 6 is even more noticeable than that 
on (/>'s period. 

In the same way, it's easy to derive A6' has unit magnitude when 7^ ~ m/^/X\(f>\ <^ 
1. Normally in an expanding universe where the backreaction is negligible, \5'\ is always 
decreasing so the resonance would gradually undergo a shift from the large 5' stage to the 
small 6' stage. The relation of 7^ however tells that the parametric resonance may be subject 
to a shift back to the large 6' stage from the small 6' stage, or even, the small 5' stage may 
not appear at all. The back shift, whenever exists, always lies ahead of the explosive boost. 
So for deep discussions of the ending conditions of the exponential boost, one should always 
take it into account. 

4 Cosmological Applications 
4.1 Preheating 

Preheating is a possible stage after inflation, during which the inflaton decays into particles 
through parametric resonance. Preheating was first realized in [15, 16], and the parameter 
space is the same with what we have discussed. Therefore, during preheating, the number 
of preheated particles grows exponentially, much faster than the redshift rate from universe 
expansion. 

Although preheating should be followed by particle thermalization, and possibly a sec- 
ond stage of reheating between them, they are not within our consideration in this paper. 
Instead of calculating the reheating temperature and how long reheating and thermalization 
last, this section is devoted to showing how this framework may be applied on specific pre- 
heating calculations. Details about the second stage of reheating and thermalization can be 
found in [16] and [22]. For the same reason, we don't care the reheating ratio between the 
dark and visible sector. To solve the dark radiation problem, one might want to refer to 
inflation models with more particle physics foundation, e.g. SUSY inspired inflation and the 
subsequent reheating and thermalization studies [23-25]. 

During preheating, the universe is dominated by (j) in our model, which is usually the 
inflaton. As a demonstration, we here just consider the simple case of the preheating of 
a generic slowroll inflation with inflaton (p. To solve problems like horizon and flatness, 
inflation requires an e- folding at least ~ 60. At the beginning of such slowroll inflations, 
there should be Ax^ <^ m?. Otherwise there would be a signiflcant decrease in c/i's effective 
mass and thus no inflation although (j) would still be slowroUing. 
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In general, the end of slowroll inflation is signaled by the destruction of the slowroll 
condition of (p, after which <j) begins to oscillate and preheating takes place. We would like 
to neglect the first cycle or so of (p in the beginning of preheating here, so ^ H holds 
for the rest of the preheating period. Because Ax^ <^ m? holds in the beginning of inflation, 
and X is redshifted greatly during inflation, this relation is expected to still hold for a very 
long period in preheating. Because the backreaction terminates parametric resonance once it 
comes to important, it acts as an ending condition of the preheating scenario. Therefore we 
can first neglect the presence of backreaction, and then calculate the take-over of backreaction 
separately. 

We label the beginning of inflation as 0, the end of inflation (i.e. the start of preheating) 
as 1, and the point we start to consider in preheating as 2. The evolution of the universe is 
then divided by such points, whose labels we use as subscripts, into intervals such as 1 — t- 2 
for the convenience of future references. To be exact, we define the point 2 as 

\H=e'^\<Pi\- (4.1) 

X is redshifted during inflation (0 — )• 1), but the boost and redshift effect approximately 
cancel during 1 — )• 2. We therefore get 

|X2| ~ Ixil = e"5^|xo|- (4.2) 
The universe is dominated by during preheating so the Hubble rate writes 

^^ = ^=F (4.3) 



3M2 

in which for point 1 and after, 

\ct>\ = \ct>i\[^f. (4.4) 

After point 2, the system starts at the large 5' stage. There may be a transition from 
large to small 8' stage during the evolution, but the termination of parametric resonance by 
backreaction is certain to take place at the large 6' stage, which is ensured by its influence 
on 5' shown in eq. (3.27), slightly before the positive feedback ends the exponential boost. 
We thus also label the point after which the large 5' is ensured by backreaction, and the 
point backreaction becomes important and consequently terminates parametric resonance as 
points 3 and 4, also in chronological order. 

After defining the timeline, we now work on the evolution of x- The large and small 
5' stages share the same average boost rate (r/), so after point 2, field x is affected by 

parametric resonance and becomes proportional to e~™^*~*^^. Due to universe expansion, 

Ixl is redshifted and thus acquires the factor (02/0)2. As we have demonstrated in section 

-i 3 
2.1, Ixl oc rn^^ , this gives from the redshift of (j) field, another factor (0/02)4. We 

multiply them all together to get the statistically average evolution of field x 

lxl = IX2|(^)^e^™(-^). (4.5) 



In a matter-like universe, we have for any time ta and tf,, 



3 



a{tb)Y 1 J. 67rp(ta) 
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which, under the condition loga(tfe) — log a(ta) > 1, can be simphfied to 

a{tb) \ ^ Mp 



where p{ta) is the total energy density of the universe at time ta- Because (p dominates 
preheating, the evolution of |x| finally writes 

lxl = IX2|(f)^^^^^^\ (4.8) 

in which 

{3tt)-2\<P2\ 

From the above equation, we immediately see |x| is actually independent of A and m. 
It's obvious \x\ is A-independent from the constant boost rate. The dependence on m, the 
bare mass of (p, is however canceled by the same m-dependence of p. This can be also viewed 
as m defines a time scale which both the Hubble rate and the boost rate of x per time interval 
are proportional to. 

We then derive the evolution of x at points 3 and 4 which are defined from backreaction. 
Here we assume (p is the inflaton of a generic slowroll inflation of N e-folds. At the beginning 
of it, i.e. point 0, there should be A|xoP < to allow inflation to happen. After being 
redshifted all along inflation, |xi| then becomes extremely small and this allows a safe period 
of preheating before the backreaction becomes important. After that, backreaction may 
get strong and, points 3 and 4 may be met in turn at any moment. To show the relation 
'^IXoP < "i^; we deploy the parameter 7q = X\xo\'^/m? < 1. Here we want to give a 
conservative and natural estimation, so we assume —100 < log 70 < 0. 

As we have demonstrated in section 3.2 that backreaction ensures the large 5' stage at 
7^ ~> m/^/X\(f>\, we therefore define point 3 as 



2 -73 = ^^- (4.10) 



Given the evolution of both fields, eq. (4.4) and eq. (4.8), this can be reduced to the equation 
of A'^2-s>3) the e-folding the universe experienced between point 2 and 3, 

gei^^^3 _ jV2^3 = AT + ^ log ^ . (4.11) 

Also, we can work on point 4 in the same way. After defining point 4 as A|x4p = ""t-^, 
the calculation is straightforward to give 

Bel^^^^ - ^iV2^4 = ^ - 3 log 70- (4.12) 

In both of the above equations, the dominant term on the l.h.s is the exponential one. 
We can thus see the number of e-folds of preheating is only very weakly dependent on A and 
70, because they are inside two "log"s. This is because (j) oscillates faster and faster as the 
universe expands, and in consequence x acquires a much larger boost rate per Hubble time 
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Figure 10. N2-^3 and N2^i as functions of A. Blue and purple correspond to A^2-i.3 and A^2^4 
respectively. It's easy to see on the figure that A'2-i.3 is only weakly dependent on A and N2^a is 
independent of A. Parameters are chosen to accord with the quadratic slowroU inflationary model. 



ends at \4>i\'^ = 27Mp/167r, and has the number of e-folds of inflation N = 60. We also need 
the value m ~ IQ-^Mp, which comes from the strength of observed CMB anisotropy. Figure 



To provide a brief example, we take the simple case of quadratic slowroll inflation which 

10 shows the relations N2^3 and N2^4, have with A at constant 70 = 0.1. The dependence 
of 7o is weak, and a smaller 70 simply shifts both A''2-5>3 and N2^a upwards a bit with the 
same amount. For 70 = e"^"'', there is N2-^4 = 4.86. 

From eq. (4.12), we can also calculate the value of the r.h.s from For N2^4 = 4 

and 5 respectively, there is — | log 70 = 32 and 155. As a result, for most cases where 
is not very large and 70 is not too small, preheating is typically ended by backreaction at 
the 5th e-fold of parametric resonance i.e. the 6th e-fold after inflation, provided no other 
mechanism terminates preheating before that. 

In section 3.2, we have suggested the correction from backreaction to 6' always comes 
to important earlier than that to m^. In figure 10, it agrees with large A. For small A, 
the figure suggests a reverse in their sequence. The reverse however actually reflects the 
violation of A|(/)p ^> m^. For values of A on the left of the intersection point of the two curves 
in figure 10, the condition X\(j)\'^ ^> should be broken well ahead of point 4 being reached, 
so the behavior of the system afterwards in fact can't be predicted by our theory. This 
indicates another possible exit of parametric resonance, which is caused by a small \(p\ that is 
comparable to m/y/X. We label this exit as subscript a, and from the relation A|(/>ap = m^, 
we can easily derive the e-folding 



iV2^a=^log^. (4.13) 



1, MH^ 

-log 

Since both point a and 4 act as exits of parametric resonance, it is important to know 
which point comes earlier to find out by who and when parametric resonance is terminated. 
By setting N2^4: = N2^a, calculations give a critical value for A 

(A^- |log7o)2m2 , 



- 22 - 



For X > Xa, parametric resonance ends with backreaction. The ampHtudes of (p and x reach 
the same scale of magnitude in the end and equihbrium may be reached between the two fields. 
For X < Xa however, parametric resonance is terminated by the violation of A|0p ^ m^. In 
this case, the evolution afterwards is not characterized by our theory, and we only know 
at the time this condition is broken, equilibrium has not yet been reached. Moreover, for 
A < m?/\4>2\'^, there is N2^a < 0. This indicates the condition of our framework is not 
satisfied and such a parametric resonance does not take place immediately after inflation. 

Meanwhile, we are also interested in the transition point from the large 5' stage to the 
small 5' stage. In the context of preheating, 5' is defined as 



When 5' reaches e ^ ~ 0.1, we think the transition takes place and label it as point h. We 
thus derive the e-folding from point 2 to 6 



For N2^b > -^2-5>ai parametric resonance is terminated before the transition point, so 
the whole process is large 5' . Similarly if A''2-s.f, > A''2-s.3, the large 5' stage is also guaranteed 
so no small 5' stage takes place. For N2^b < 0, parametric resonance starts from the small 5' 
stage, and whether a large 5' stage presents (although very short) depends on how preheating 
is terminated. 

To conclude this subsection, we perform a more precise calculation of the parametric 
resonance of quadratic slowroll inflation, and give figure 11 which divides the parameter 
space of A and 70 into five sections. The small 5' stage may either exist or not exist for 
the terminations by both backreaction and the small \(f)\. There is also the case that no 
parametric resonance shall take place within our framework, so in all there are five sections. 
The numbers of the total e-folds of parametric resonance and the e-folds of the small 5' stage 
are also contoured in figure 11. For the termination by small \(j)\ however, the reheating 
is incomplete and there should be a subsequent decay of inflaton after the termination of 
parametric resonance, but our theory is unable to characterize it. From the figure, we also 
confirm that if preheating is terminated by backreaction, it should typically happen at the 
5th e-fold of preheating, i.e. the 6th e-fold after infiation. 

4.2 Locked Infiation 

Locked infiation, or new old infiation[9], is a fastroU inflation scenario during which a fastroll 
field (j) locks the other field x ifs false vacuum. The false vacuum has a nonzero energy 
density which drives infiation. Locked inflation has the same potential shape with hybrid 
inflation [5], but the parameter configuration is different. 

Not long after the proposal of locked inflation, however, opposite opinions are posted 
suggesting underlying problems[20, 26]. In [20], they concluded that saddle inflation, loop 
correction and parametric resonance problems eliminates part of the parameter space each, 
while they three together kill all possible parameter conflgurations. At least one of those 
problems therefore has to be fixed to make locked infiation work again. During locked 
infiation, parametric resonance is disfavored for its exponentially growing energy transfer 
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Figure 11. The different circumstances of prelieating in tlie parameter space of A and 70. The 
green and red regions are backreaction termination with and without the small S' stage respectively. 
The blue and orange regions are terminated by the small |(/)|, with and without the small S' stage 
correspondingly. In the gray region, the preheating stage is not entered after inflation. Among the 
dashed lines, the number of total e-folds of parametric resonance is contoured in black and that of 
the small S' stage is contoured in red. 

from (j) field to X) so it acts as a problem here which may end locked inflation much earlier 
than expected. 

Prom our above analysis of the mechanism of parametric resonance, we infer for (j) with 
more than one real component, parametric resonance is inhibited. This solves the parametric 
resonance problem of locked inflation and makes it possible. Such an argument comes from 
eq. (2.5), the condition of rolling stage. For (j) singlet, is dominated by Xcj)'^ so this 
condition breaks when (p is approaching zero. Only when eq. (2.5) breaks, (p would enter the 
zero-crossing stage and perform the exponential boost. If Xcp'^ isn't the absolute dominant of 
m^, eq. (2.5) would never break down and the zero-crossing stage would never be entered, 
even if (/> = 0. For (p multiplets, every component of (/> interacts with x independently. So 
when one component crosses zero, other components are in general nonzero and they provide 
X with a large effective mass. The zero-crossing component thus never dominates and eq. 
(2.5) never breaks. 

As an example, we take a (p doublet, whose components are (pi, (p2- For demonstration, 
we simply assume they have an identical amplitude \(p\ but a phase difference | — 2(/?. So 
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they can be written as 

= \(j)\cos{m^t - ip), (4.17) 
(p2 = \4>\sin{m^t + ip). (4.18) 

The effective mass squared of x then becomes 

= M^ + 2A(0? + (1)1) = + 2A|(/>p(l + sin2(/jsin2m^t), (4.19) 

equivalent to the singlet case with transformed — )■ M'^ = + 2A|(/)p, A — )• A = A sin 2^9 
and rrifj, — )• = 2m0. 

The equivalent singlet case however has the relation X\(j)\'^ < M^, which breaks the 
assumption A|(/)p » » M^. In such cases, the condition of rolling stage eq. (2.5) would 
never break, and the parametric resonance effect would then be inhibited. 

In general cases where \(pi\ ^ \(p2\ but \(pi\ ~ \<p2\, there are also similar results, which 
typically have ^ 0.1A(|(/)i|2 + \(p2\^). For the same reason, eq. (2.5) would never break 
and locked inflation can last long enough for cp multiplets. Therefore for (p multiplets, the 
parametric resonance "problem" of locked inflation is not actually a problem. Of course, 
this doesn't solve the problem of embedding locked inflation in some particle physics theory. 
However, one still needs a viable particle physics setup (or a convincing argument about the 
decay to visible sector) to make locked inflation promising. This is beyond the coverage of 
this paper. 

5 Discussion of Inhomogeneity Amplification 

The cosmic inhomogeneity, if arose from the inflaton, may get amplified by parametric res- 
onance during the large 6' stage. The amplification also requires a change in the equation 
of state of the universe that is triggered directly or determined indirectly by the total boost 
rate of \x\- This process is similar to modulated perturbations discuss in [27]. The difference 
here is the decay rate is determined by the relative phase difference between the infiaton and 
the preheated field, i.e. both fields instead of the infiaton only. 

To explicate the whole process, we start with two points at a distance 1/k. In the 
beginning of inflation, they live within a Hubble radius, so homogenous x is assumed. After 
the mode k leaves horizon, the two points become uncorrelated and the phases of x at the two 
points start to differ because of the inhomogeneity 5(j)k which causes different ra^^s. Also, as 
we choose the end of inflation at a speciflc value of 4> on the uniform-density hypersurface, the 
point with a slightly larger would also gain a slightly longer inflation. For both reasons, at 
the point where (j) is larger, the phase that x experiences would also be larger. This transfers 
inhomogeneity to x's phase from energy density. 

After inflation, the universe enters the large 6' stage. Due to statistical reasons, the 
total transformation matrix of the whole large 5' stage would be an ellipse with the aver- 
age total boost rate (77) and the difference between the major and minor axes of the scale 
2A(?]). Different initial phases of x then correspond to different total boost rates, i.e. dif- 
ferent parametric resonance efficiencies. If there exists a transition of the equation of state 
depending on the efficiency of parametric resonance, like the backreaction end which changes 
the universe from matter-like to radiation-like when A|xP ~ w-^, additional inhomogeneity 
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may be generated."' At points where the efficiency is lower for example, more time is needed 
to boost X to the amplitude required for the transition in the equation of state to take place, 
and subsequently the entrance into a different equation of state is delayed. Therefore there 
are regions which stay in the larger equation of state longer and suffer a heavier redshift 
due to universe expansion, and also regions which enjoy a weaker redshift. The difference in 
redshift rates finally leads to different energy densities at different places, and thus additional 
curvature perturbation on the uniform-density hypersurface. If the effect is stronger than 
the existing inhomogeneity, we will find inhomogeneity amplified in this process. 

To calculate, we start from inflation in which the inhomogeneity transfers from (j) to 
the phase of x- We still adopt the two points at the distance 1 /k and difference in (p at the 
horizon exit of mode k as (50^. Because we use a universal ending condition for inflation, we 
align the evolution history of the two points at the end of inflation. Then the only difference 
5^k causes to the phase of x is one point enjoys an additional period of inflation as the 
field (j) rolls from cj) + 5(i)k to and this gives an additional phase to x- Using the slowroll 
approximation, the resulting phase difference can be calculated 



A99 = V2A 



/3A 3Afpff| 5^ 



where we use the subscript k to represent values at the horizon exit of mode k. 

During the large 6' stage, the parametric resonance process is approximated as stochas- 
tic. So we acquire the statistical result of the average total boost rate (77) and its statistical 
error A{ri), which are defined in eq. (3.12) and eq. (3.14) and are functions of n, the number 
of half cycles of (j) during the large 5' stage. If we consider the total transformation matrix 
of the whole large 5' stage which transforms a circle of unity to an ellipse, (77) can also be 
understood as the typical average total boost rate of the ellipse and A (77) as half of the typical 
difference of the total boost rate between the major and minor axes. 

Now we have the ellipse as the total boost rate function of the initial phase. The total 
boost rate goes from maximum to minimum with a change of initial phase tt/2, so the phase 
difference Aip would give a difference in the total boost rate, and thus a difference in the 
number of half cycles of (p needed to reach the same amplitude For A(p < 1, we have the 
number of additional half cycles 

2Aw 2A(?7) , , 

An=—^—^. 5.2 

Then An leads to delayed or advanced transition of the equation of state. Let us suppose 
the equation of state is changed by Aw. The resulting inhomogeneity is then 

6pk_An7r [6X ISAwMpHHj Ajrj) 6<Pu 

Pk m V vr 771'' (77) (pk 

where H here is the Hubble rate at the transition of the equation of state. 



^During preheating, the universe is matter-hke because 4> dominates with its energy density proportional 
to . For preheating ended by backreaction, the dominant potential term becomes A^^x'^- Therefore (j) and 
X have the same redshift rate. With their redshift rate for constant-mass fields, and their scaling relation 
with effective mass eq. (2.11), we can conclude both fields are proportional to . This gives the energy 
density A|(^p|xl^ oc and thus a radiation-like universe. This result is also numerically verified in [28]. 
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If we don't want it to overwhelm the existing inhomogeneity 5(f)k/(t>k^ we then need it 
to be smaller than that, which gives 



and consequently, 



For > 1, we have instead 
which generates the inhomogeneity 



An=?^, (5.6) 



6pk GttAwH A{ri) 
Pk m {i]) 



(5.7) 



A{'q) / {rf) is at least 1 and increases as the large 6' stage continues; Aw ~ |- To make 
^Pkl Pk < 10^^, we therefore need at least 8 e-folds before the transition to get a small enough 
-ff , which is hardly possible for a quadratic slowroll inflation, unless we have thousands of 
e-folds of inflation to redshift x- 

So we will still focus on the Aip < 1 case. If we consider the cosmic scale whose 
^i'k/'t'k ~ 10~^ and ~ lOm^, we can see from eq. (5.1) the constraint from Aip < 1 
is already quite strong. For the backreaction end of preheating of the quadratic slowroll 
inflation, there is A < 10~^. 

However the constraint from eq. (5.5) is even stronger. At the backreaction end, the 
equation of state switches from to | , so Aw = | . For the preheating fully consisted of the 
large 6' stage, we can simplify part of eq. (5.5) to a function of A'^2-!>4 only 



{r]Y 3V37rMp(r?) 



e^^^-*. (5.^ 



m A{r,)^ 4|<A2|((??2)-(7?)^ 



So eq. (5.5) becomes 



2592Mp|,/.2|i/,Mr/2)-(r/)2' ' ^^'^^ 

Here we simply choose N2~^i = 4.5 and H"^ = lOm^ to get the numerical value A < 10~^^. 
If we want the result A < 10~^, from which we can reserve a tiny living space for the 
backreaction end as shown in figure 11, we would need A^2-s>4 > 18, which means more than 
10^^ e-folds of inflation. For the backreaction end with both the large and small 5' stages 
during preheating but the same e-folding of the large 5' stage, the constraint is even stronger. 

Such a strong constraint basically rules out the backreaction end of preheating for the 
quadratic slowroll inflation, unless additional techniques are taken. The attempt of using a 
dynamic A to solve this problem is also difficult, because that would require A < 10"^^ during 
inflation, which makes x ^-Iso slowrolling. For natural initial conditions, one then needs to 
deal with a very large x held after inflation which doesn't even need parametric resonance 
to become large. This however disagrees with our initial purpose of choosing x ^ the first 
level decay product of cj). Therefore, we infer that the first level of decay of the inflaton of 
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quadratic slowroll inflation is insufficient for reacliing a cliemical equilibrium between the two 
types of particles, if we still adopt 5(1) /(I) ~ 10~^ at the cosmic scale. 

One may, however, assume inflation only generated a much weaker inhomogeneity and 
preheating then amplified it to the proper amount. To get a weaker inhomogeneity, we need 
to change the mass m. In order to benefit from the scaling property, we will still use mo 
to indicate the old m which alone can get the proper inhomogeneity during inflation. On 
the r.h.s of eq. (5.3), H and are both proportional to m, and S(j)k/4>k oc m^, so the r.h.s 
as a whole is proportional to m, while the l.h.s is the desired inhomogeneity and should be 
kept unchanged. We can absorb the dependence on m into A for the rest of the calculation, 
and drag it out in eq. (5.5). If nothing else is modified, we would get a similar result 
Am^/mQ ~ 10""'^'^. Still if we want A ~ 10~^, this would require m ~ lO^-'^^Afp. Also, because 
the amplification rate is the term on the l.h.s of eq. (5.3) before 6(pk/4'k which has in it, 
the amplification is scale dependent. Different length scales leave horizon at different Hubble 
rates which consequently give different inhomogeneity transfer rates. This comes from the 
different amounts of phase differences of X; which are only generated at their corresponding 
horizon exits. Therefore, the scale dependence further adds — 2e to the spectral index for this 
case, providing e = g^-^ is the slowroll parameter. 

For generic cases, we just need to redefine Aip. Still aligning every point at the end of 
inflation and assuming the only difference the inhomogeneity of inflaton gives on x's phase 
is a slightly longer or shorter inflation, we can then define it as Aip = {my.)At. Here {rriy.) 
is the average during the additional inflationary period At caused by the inhomogeneity 
of the inflaton. For short, we use 6k = Spk/p as the energy density inhomogeneity at horizon 
exit. The time difference of inflation can then be written as 

dp dlogp 2Hkek' 

where again the subscript k indicates the value at horizon exit of mode k. So the inhomo- 
geneity arose from parametric resonance, represented as 6k, should be 

6k = T\ — u ^f"- ^-^^ 

The amplification rate then should be the terms on the r.h.s before 5k (for Aip < 1). 

Although we have only given one example, this constraint has basically ruled out the 
backreaction end of parametric resonance in preheating for m ~ 10~^. This constraint is 
also quite strong in general cases, as long as all of the following criteria are met so that the 
constraint is applicable. 

1. The inflaton has some inhomogeneity after the horizon exit. 

2. There is a period of large 6' stage during the flrst level decay of the inflaton. 

3. A transition of the equation of state of the universe is a direct or indirect consequence 
of the flrst level decay product's reaching a critical density. 

This constraint however doesn't apply to the small 6' stage, because all phases converge 
to the same one in the small 5' stage. Nor does it apply to perturbative reheating, etc. The 
constraint doesn't affect models without any transition in the equation of state, such as the 
inflaton behaves like radiation after inflation which has the same equation of state with the 
interaction term. 
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6 Summciry 



In this article, we have estabhshed an analytic framework for the homogeneous mode of 
parametric resonance in a particular parameter space. In a static background, the Mathicu 
equation is analytically solved to give the stable and unstable bands. In an expanding 
universe, distinct stages of large and small 5' are separately considered. At the large 5' stage, 
the phase x experiences for every half cycle of (f) is regarded as random and independent 
of the history, like the Markov chain. During the small 5' stage, all states with different 
phases converge to the same eigenstate which is slowly varying. In this eigenstate, |x| enjoys 
a stair-like growth. Both stages share the same average boost rate (77) = 0.346. The effect of 
backreaction is considered, in which a positive feedback acts as the terminator of parametric 
resonance once it takes over m^. We also find the backreaction has an impact on 5' which 
comes to important earlier than the take-over and guarantees a short large 5' stage before 
the backreaction exit of parametric resonance. 

Applying it on preheating, we have constructed a timeline starting from the beginning 
of inflation till the backreaction exit. Taking the quadratic slowroll model as an example, 
the large 5' guarantee point and the transition point from large to small 5' are included. The 
violation of A|^p ^ w? has also been considered as another exit of parametric resonance, 
namely the small 101 exit. Regarding to how parametric resonance is terminated and whether 
a small 5' stage presents, we divide the parameter space of A and 70 into four parts. The 
backreaction exit typically takes place during the 5th e-fold of preheating for slowroll models 
with e-folding 60. 

However in a subsequent discussion, we have found parametric resonance during pre- 
heating can be responsible in amplifying inhomogeneity. To avoid conflict with current CMB 
spectrum, it acts as both a constraint on parameter space for some inflationary models, and 
an amplifier to the weak inhomogeneity of some other models which are otherwise disfavored. 
We have also set up a demonstration by ruling out the backreaction end of preheating of the 
quadratic slowroll inflation with mass 10~^. Amplification rate for general cases is calculated, 
and criteria of applicability of this effect is explained. 

On the other hand, the parametric resonance problem of locked inflation isn't a problem 
for (p multiplets. This is because when one component of ^ crosses zero, the contribution 
from other components on would be strong and consequently prevent the exponential 
boost effect. 

This paper aims to show an alternative framework for the broad resonance. For this 
reason, we don't specify the particle physics model and only use singlets. We also show the 
preheating stage for quadratic slowroll inflation only as an demonstration, without calculating 
the reheating temperature and thermalization, or discussing difficulties of inflationary models, 
such as the visible v.s. dark sector or how to embed locked inflation in a particle physics 
foundation. 
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